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w  IS.  ABSTRACT 

P  T-his  report  describes  the  results  of  a  theoretical  study  on  the  analysis  and  design 
of  leading  edge  slat  plus  main  airfoil  combinations.  Originally  it  was  planned 
fee  approach  this  problem  using  three  basic  methods:  X--Jtepresentation  of  the 
slat  by  a  point  vortex,  if.  Representation  of  the  slat  by  a  set  of  thin  airfoil 
theory  singularity  distributions,' 4ind  III.  exact  solution  of  the  complete  two 
element  airfoil  problem  using  conformal  mapping  methods.  After  developing  the  basic 
analytical  work  for  method  III,  it  became  apparent  that  the  remaining  analysis  and 
computer  programming  required  for  even  a  rudimentary  evaluation  of  this  approach 
wpuld  exceed  the  total  time  available  for  the  entire  study.  Therefore,  it  was 
''decided  to  concentrate  on  the  development  of  Methods  I  and  II.  An  "on-line"* 
(Computer  ggapLj  c.?  progranrwm  prepared,  for  solution  and  evaluation  of  Method  I* 
JThe  main  airfoil  l*  chosen  from  an  infinite  set  of  airfoils  whose  mapping  functions 
to  the  unit  circle  are  analytically  defined.  (The  classical  Joukowski  airfoils 
are  members  of  this  set.)  The  vortex  strength  and  position  with  respect  to  the 
main  airfoil  is  then  prescribed  and  tue  program  calculates  and  displays  the  velocity 
distribution  on  the  airfoil  with  ana  without  the  vortex  present.  This  program  lias 
I  been  applied  to  the  study  and  development  of  tome  basic  slat  design  guidelines 
which  are  presented.  The  problem  of  designing  A  leading  edge  slat  which  provido6 
a  specified  velocity  modulation  on  the  nose  region  of  a  main  airfoil  hao  baan 
studied  using  the  approach  of  $kfthod  II.  ^A  computer  program  has  been  prepared 
which  for  a  given  slat  thickness  distribution  designs  a  slat  camber  line  that 
!  provides  a  specified  velocity  modulation  on  ti\e  nose  region  of  an  ellipse.  Sample 
results  are  included  for  a  limited  set  of  test  cases. _ 
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1.0  ABSTRACT 

This  report  describes  the  results  of  a  theoretical  study  on  the  analysis  and 
design  of  leading  edge  slat  plus  main  airfoil  combinations.  Originally  It 
was  planned  to  approach  this  problem  using  three  basic  methods:  I.  Repre¬ 
sentation  of  the  slat  by  a  point  vortex,  II.  Representation  of  the  slat  by 
a  set  of  thin  airfoil  theory  singularity  distributions,  and  III.  Exact  solu¬ 
tion  of  the  complete  two  element  airfoil  problem  using  conformal  mapping 
methods.  After  developing  the  basic  analytical  work  for  Method  III,  It  be¬ 
came  apparent  that  the  remaining  analysis  and  computer  programming  required 
for  even  a  rudimentary  evaluation  of  this  approach  would  exceed  the  total 
time  available  for  the  entire  study.  Therefore,  It  was  decided  to  concentrate 
on  the  development  of  Methods  I  and  II. 

An  "on- Vine"  computer  graphics  program  was  prepared  for  the  solution  and 
evaluation  of  Method  I.  The  main  airfoil  Is  chosen  from  an  Infinite  set  of 
airfoils  whose  mapping  functions  to  the  unit  circle  are  analytically  defined. 
(The  classical  Joukowskl  airfoils  are  members  of  this  set.)  The  vortex 
strength  and  position  with  respect  to  the  main  airfoil  Is  then  prescribed  and 
the  program  calculates  and  displays  the  velocity  distribution  on  the  airfoil 
with  and  without  the  vortex  present.  This  program  has  been  applied  to  the 
study  and  development  of  some  basic  slat  design  guidelines  which  are  pre¬ 
sented  in  this  report. 

The  problem  of  designing  a  leading  edge  slat  which  provides  a  specified  velo¬ 
city  modulation  on  the  nose  region  of  a  main  airfoil  has  been  studied  using 
the  approach  of  Method  II.  A  computer  program  has  been  prepared  which  for 
a  given  slat  thickness  distribution  designs  a  slat  camber  line  that  provides 
a  specified  velocity  modulation  on  the  nose  region  of  an  ellipse.  Sample 
results  are  Included  for  a  limited  set  of  test  cases. 
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5.0  INTRODUCTION 

The  designer  of  high  lift  systems  for  modem  aircraft  Is  faced  with  a  severely 
over-constrained  problem.  The  satisfaction  of  cruise  configuration  geometry 
together  with  structural  and  mechanical  requirements  quite  often  results  In 
a  high  lift  system  which  Is  less  than  optimum  aerodynamlcally.  An  additional 
complication  results  from  the  fact  that  even  If  the  structural  and  mechanical 
requirements  were  relaxed,  the  aerodynamic  design  of  high  lift  systems  Is 
still  quite  empirical  and  heavily  reliant  on  wind  tunnel  testing.  It  Is  the 
purpose  of  this  study  to  attempt  to  Increase  the  knowledge  and  understanding 
of  one  of  the  basic  components  of  a  high  lift  system,  the  leading  edge  slat. 

A  primary  function  of  a  leading  edge  slat  Is  to  suppress  the  leading  edge 
pressure  peak  on  the  main  airfoil  when  It  Is  operated  at  high  angles  of  at¬ 
tack,  and  thus  avoid  boundary  separation  and  possible  stalling  of  the  main 
airfoil.  Referring  to  figure  1,  the  circulation  about  the  slat  creates  a 
counter  velocity  component  In  the  nose  region  of  the  main  airfoil  which  re¬ 
sults  In  a  reduction  of  the  pressure  peak. 

Although  simple  In  principle,  the  design  of  efficient  and  effective  leading 
edge  slats  Is  a  very  complex  problem.  In  fact.  It  Is  not  particularly  clear 
that  a  "best"  slat  design  exists  or  can  even  be  rationally  defined.  On  the 
other  hand,  a  slat  which  provides  a  desired  Cj_  capability  while  satisfy¬ 
ing  the  mechanical  and  structural  requirements  Is  certainly  a  successful 
design  solution.  In  this  light,  some  characteristic  aerodynamic  slat  design 
problems  can  be  defined  whose  solutions  should  provide  a  significant  advance 
In  the  understanding  and  design  capability  of  leading  edge  slats.  It  was 
not  Intended  that  all  of  the  following  problems  be  approached  and  solved 
under  the  limited  scope  of  the  present  study,  however,  they  are  briefly  dis¬ 
cussed  here  In  order  to  provide  perspective  for  the  present  work. 

(a)  Design  of  a  slat  which  separates  from  the  leading  edge  of  the  main 

airfoil. 

This  Is  probably  the  most  common  of  the  slat  design  problems,  and  also 
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the  most  difficult  from  the  aerodynaml  cist's  point  of  view.  Here  the 
main  airfoil  geometry  Is  already  specified  and  the  designer  Is  asked  to 
separate  a  portion  of  the  leading  edge  to  form  a  slat  as  shown  In  figure 
2.  The  upper  surface  geometry  of  the  slat  Is  fixed,  and  the  lower  sur¬ 
face  of  the  slat  and  the  slat  and  the  nose  of  the  main  airfoil  must  be 
shaped  so  they  will  fit  together  when  the  slat  Is  retracted  as  shown  In 
figure  2.  The  resulting  slat  Is  then  oriented  with  respect  to  the  main 
airfoil  to  provide  the  desired  performance. 

(b)  Design  of  a  slat  which  provides  a  specified  modulated  pressure  distri¬ 
bution  on  the  nose  region  of  a  given  main  airfoil. 

This  problem  Is  less  constrained  than  problem  (a)  above.  The  slat  geo¬ 
metry  and  orientation  are  left  completely  free  and  It  Is  only  required 
that  the  resulting  slat  provides  the  specified  pressure  modulation  dis¬ 
tribution  on  the  nose  of  the  main  airfoil.  In  addition,  the  flow  should 
not  separate  on  the  slat  Itself. 

(c)  Simultaneous  design  of  both  the  slat  and  the  main  airfoil. 

It  Is  obvious  that  this  problem  requires  at  least  some  generalized  con¬ 
straints  or  performance  goals  to  have  any  solution  at  all.  However  a 
capability  In  this  area  could  provide  solutions  to  more  specific  problems 
such  as:  "design  a  two  element  airfoil  system  which  obtains  the  maximum 
possible  lift  without  separation". 

Problems  (a),  (b),  and  (c)  above  are  of  course  not  distinct,  nor  do  they 
represent  In  themselves  the  entire  spectrum  of  leading  edge  slat  design  prob¬ 
lems.  They  do  serve  to  Illustrate  the  basic  forms  which  such  problems  take 
on,  and  this  suggests  the  types  of  analyses  which  are  most  suitable  In  attempt 
Ing  to  obtain  design  solutions. 

There  exist  several  theoretical  analysis  methods  which  may  be  applied  to  the 
leading  edge  slat  design  and  analysis  problems  just  described.  One  of  the 
most  powerful  and  versatile  methods  for  the  analysis  of  a  multi -element 
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airfoil  system  of  specified  geometry  Is  the  Douglas  Neumann  potential  flow 
computer  program  (reference  1).  This  program  provides  exact  solutions  for 
the  potential  flow  field  (on  and  off  the  bodies)  for  multi-element  airfoil 
systems  of  up  to  five  elements,  and  therefore  Is  an  extremely  valuable  tool 
for  evaluating  and  predicting  the  aerodynamic  performance  of  leading  edge 
slat  designs.  This  type  of  solution  Is  called  "direct"  In  that  It  requires 
the  airfoil  and  slat  geometry  as  an  Input.  It  can  be  used  as  a  design  tool 
by  examining  the  results  for  a  particular  geometry,  and  then  adjusting  the 
geometry  and  re-running  the  program  until  the  desired  performance  is  obtained. 
Unfortunately  the  determination  of  what  geometry  adjustments  should  be  made 
requires  considerable  experience  and  Intuition,  and  this  procedure  tends 
to  be  rather  costly  and  time-consuming.  However,  this  approach  appears  to  be 
the  best  available  for  problem  (a)  above  at  this  time. 

An  analytic  procedure  for  Iterating  with  direct  solution  results  has  been 
developed  by  Wilkinson  in  reference  2  for  the  design  of  multi-element  air¬ 
foil  systems.  The  pressure  distributions  on  the  upper  surfaces  of  the  air¬ 
foil  elements  are  prescribed,  and  an  Initial  geometry  Is  specified.  Wilkinson's 
Iterative  procedure  then  systematically  modifies  the  camber  distributions  of 
the  airfoil  elements.  The  maximum  thicknesses,  thickness  distributions,  chord 
lengths,  and  gaps  between  airfoil  elements  are  maintained  at  all  stages  as 
equal  to  the  respective  quantities  on  the  Initial  geometry.  The  advantage  of 
this  type  of  method  lies  In  the  fact  that  If  a  final  solution  Is  obtained 
(i.e.  It  converges)  It  Is  exact.  However,  there  is  no  guarantee  that  an 
arbitrarily  Input  set  of  Initial  geometries  and  upper  surface  pressure  dis¬ 
tributions  can  be  iterated  to  convergence.  This  method  Is  useful  for  approach¬ 
ing  the  solution  of  problem  (c)  above. 

O'Pray  In  reference  3  has  developed  a  linearized  solution  for  problem  (b)  which 
for  a  given  slat  chord,  slat  orientation,  and  slat  thickness  distribution, 
shapes  the  slat  camber  line  to  provide  the  specified  pressure  modulation  dis¬ 
tribution  on  the  nose  region  of  the  main  airfoil.  It  appears.  In  principle, 
that  this  method  can  be  extended  to  obtain  a  partial  solution  to  problem  (a): 
namely,  Instead  of  specifying  the  slat  thickness  distribution,  the  slat  upper 
surface  shape  could  be  specified. 


L 


9 


MDC  J5195 


In  the  present  work,  It  was  originally  decided  to  study  the  problem  of  lead¬ 
ing  edge  slat  analysis  and  design  using  three  basic  approaches: 

I.  Simulation  of  the  slat  by  a  point  vortex. 

II.  Simulation  of  the  slat  by  a  finite  set  of  singularity  distributions. 

III.  Exact  solutions  using  conformal  mapping  methods. 

The  choice  of  these  approaches  was  based  on  the  following  reasons.  The  point 
vortex  simulation  of  a  slat  represents  the  simplest  possible  theoretical  model, 
and  Its  formulation  Is  straightforward  analytically.  The  resulting  analysis 
can  be  used  to  demonstrate  the  basic  effects  of  slat  position  and  slat  lift 
coefficient,  and  these  results  are  quite  useful  In  both  developing  a  more 
thorough  understanding  of  slat  performance,  and  in  guiding  the  development 
of  the  more  complex  studies  of  approaches  II  and  III. 

The  representation  of  the  slat  by  a  set  of  distributed  singularities  In 
approach  II  may  be  considered  as  the  next  step  up  In  complexity  of  the  theo¬ 
retical  model.  Here, thin  airfoil  theory  is  used  to  describe  the  slat  and 
the  flow  field  which  it  Imposes  on  the  main  airfoil.  The  effects  of  slat 
chord  length,  camber,  and  thickness  as  well  as  position  and  slat  lift  coef¬ 
ficient  can  be  studied  Independently  since  thin  airfoil  theory  allows  super¬ 
position.  In  principle,  this  should  provide  linearized  solutions  of  problems 
(a)  and  (b)  above. 

Approach  III  Is  by  far  the  most  elegant  and  general  of  the  theoretical  models 
considered  In  this  study.  Since  conformal  mapping  techniques  are  to  be  em¬ 
ployed,  the  final  analysis  should  provide  solutions  to  both  direct  and  Inverse 
(given  the  pressure  distribution,  find  the  corresponding  geometry)  problems. 
Once  developed,  exact  solutions  to  problems  (a),(b),  and(c)should  be  obtainable 
using  this  method. 

The  Investigations  Into  the  exact  solution  of  slat  airfoil  combinations,  by 
mapping  Into  an  annular  domain  and  the  generation  of  special  solutions  by 
mapping  two  disjoint  circles  have  been  carried  out  to  the  point  where  It  Is 
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clear  that  both  approaches  are  promising  analytically  and  appear  practical. 
However,  the  funding  under  this  contract  was  not  sufficient  for  further  com¬ 
puter  work  to  be  undertaken  to  test  these  concepts  fully,  so  that  It  was  con¬ 
sidered  advisable  to  concentrate  the  limited  effort  on  the  approximate  methods 
(I  and  II).  Nevertheless  It  Is  strongly  recommended  that  the  exact  techniques 
should  be  examined  further. 

The  following  sections  describe  the  development  of  the  analysis  of  approaches 
I  and  II.  A  limited  set  of  sample  problems  have  been  studied  and  the  results 
of  these  are  presented  and  discussed. 
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6.0  POINT  VORTEX  MODEL 


6.1  General  Discussion 


As  discussed  earlier,  the  representation  of  the  slat  by  a  point  vortex  Is  the 
simplest  possible  theoretical  model  for  studying  the  leading  edge  slat  plus 
airfoil  problem.  From  an  analytic  standpoint,  this  model  has  a  particular 
appeal  since  the  complex  potential  function  for  the  flow  field  about  a  circular 
cylinder  with  a  point  vortex  located  outside  the  cylinder  can  be  expressed 
In  a  simple  closed  analytic  form  by  using  Mllne-Thompson's  circle  theorem 
(reference  5). 

The  basic  problem  Is  described  In  figures  3  and  4  where  figure  3  shows  the 
flow  about  a  circular  cylinder  with  circulation  Iq.  Figure  4  shows  the  flow 
about  the  same  cylinder  with  the  addition  of  an  external  point  vortex  of 
strength  rs.  Merely  combining  the  complex  potential  functions  for  the  flow 
about  a  cylinder  and  for  a  point  vortex  will  not  yield  the  proper  solution. 
Instead,  an  Image  vortex  of  strength  -rs  must  be  placed  Inside  the  cylinder 
where  Its  location  Is  readily  determined  using  the  circle  theorem.  The  Image 
vortex  changes  the  circulation  about  the  cylinder  which  In  turn  moves  the 
stagnation  points  on  the  cylinder.  Addition  of  a  second  vortex  at  the 
center  of  the  cylinder  will  return  one  (but  not  both)  of  the  stagnation  points 
to  its  original  position  when  the  strength  Is  properly  set.  For  this 
study,  the  location  of  the  rear  stagnation  point  is  held  fixed  so  that  when 
the  cylinder  Is  mapped  Into  an  airfoil  the  Kutta  condition  Is  maintained. 
Therefore  the  complex  potential  for  the  flow  about  a  cylinder  with  circulation 
plus  an  external  vortex  present  can  be  written  down  In  a  simple  closed  ana¬ 
lytic  form.  It  consists  of  the  sum  of  the  complex  potentials  for  the  cylinder 
plus  those  of  the  external  vortex.  Its  Image,  and  the  circulation  correcting 
vortex. 

In  addition  to  the  classical  Joukowskl  airfoil,  there  exists  an  Infinite  set 
of  airfoil  shapes  whose  mapping  derivatives  to  a  circle  may  also  be  described 
analytically.  These  airfoils  are  referred  to  as  "pole  airfoils"  on  the  basis 
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of  the  singularities  used  to  define  their  mapping  derivatives.  A  particular 
airfoil  Is  defined  by  selecting  pole  locations,  orders,  and  strengths  within 
the  unit  circle,  e.g.  the  Joukowskl  airfoil  Is  defined  by  a  first-order  pole 
plus  a  second-order  pole  located  colncldently  inside  the  circle.  Altering 
their  location  and  strength  changes  the  camber  and  thickness  of  the  resulting 
airfoil.  A  rigorous  derivation  and  analysis  of  the  pole  airfoil  theory  to¬ 
gether  with  the  addition  of  the  external  point  vortex  Is  given  In  reference  4, 
and  therefore  only  the  application  of  this  theory  will  be  discussed  here. 

Since  the  entire  flow  field  about  a  pole  airfoil  plus  point  vortex  can  be 
described  In  closed  analytic  form,  the  study  of  this  model  lends  Itself  parti¬ 
cularly  well  to  computer  graphics.  In  this  light,  a  computer  program  using 
the  IBM  System  360  computer  together  with  the  IBM  2250  display  unit  was  writ¬ 
ten  which  provides  "on-line"  solutions.  Operation  of  the  program  Is  basically 
as  follows: 

(1)  The  strengths  and  locations  of  the  poles  which  are  to  define  the  airfoil 
are  specified  In  the  circle  plane  as  shown  in  figure  5.  (The  choice  of 
the  particular  values  used  to  obtain  a  desired  airfoil  is  described  in 
reference  4).  Figure  5  Is  a  replica  of  the  actual  display  which  appears 
on  the  IBM  2250  screen.  The  pole  specification  values  are  typed  In  on  a 
keyboard.  The  Initial  vortex  location  and  strength  Is  also  typed  in. 

(2)  Next,  a  transfer  Is  made  to  the  airfoil  plane  shown  In  figure  6  (which 
Is  also  a  picture  of  the  actual  IBM  2250  display).  The  top  figure  shows 
the  airfoil  geometry  and  location  of  the  point  vortex.  The  lower  figure 
gives  the  airfoil  surface  velocity  distribution  as  a  function  of  s,  the 
arc  length  along  the  airfoil  surface  beginning  at  the  lower  surface  trail¬ 
ing  edge  and  proceeding  clockwise  around  the  airfoil  to  the  upper  surface 
trailing  edge.  Three  velocity  distributions  are  shown:  that  for  the  air¬ 
foil  alone,  that  Imposed  on  the  airfoil  by  the  point  vortey.  <*nd  the  com¬ 
bined  result  which  represents  the  flow  about  the  airfoil  In  uie  presence 
of  the  vortex.  All  of  these  solutions  are  exact. 
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(3)  The  airfoil  angle  of  attack  and  the  vortex  strength  and  location  can  be 
changed  by  merely  keying  In  different  values  as  desired.  The  results  are 
Immediately  displayed  on  the  IBM  2250  screen  In  the  form  of  figure  6.  At 
any  time  the  display  imy  be  printed  for  later  reference  which  Is  the 
source  of  figures  5  and  6. 

(4)  The  aerodynamic  coefficients  and  parameters  listed  In  figure  6  are  defined 
and  related  as  follows: 

ALPHA  (a):  Airfoil  angle  of  attack  as  measured  from  the  chord  line. 

CL  (Cj_) :  Lift  coefficient  of  the  airfoil  alone  at  that  angle  of  attack. 

CL  TOTL  (Clj):  Lift  coefficient  of  the  total  system  (airfoil  plus  vortex). 

CL  AIRF  (Cla) :  Lift  coefficient  of  the  airfoil  In  the  presence  of  the 
vortex. 

DEL  CL  (ACl):  Specified  Increment  of  lift  for  the  total  system  (airfoil 
plus  vortex)  over  the  lift  of  the  airfoil  alone,  l.e. 

cLj  -  cL  +  ACl  (1) 

CD  AIRF  (Cq^):  Drag  coefficient  of  the  airfoil  In  the  presence  of  the 
vortex.  (Drag  on  the  vortex  Is  equal  and  opposite  to  that  of  the 
airfoil.) 

RADIUS,  THETAV  (ry,  fly):  Location  of  the  vortex  In  the  circle  plane. 

An  effective  "slat"  lift  coefficient  Is  defined  by 


-C, 


(2) 


All  of  the  aerodynamic  coefficients  are  based  on  the  airfoil  chord  and  the 
free  stream  velocity,  V^.  This  Is  more  or  less  standard  In  multi¬ 
element  airfoil  work,  and  allows  convenient  relations  such  as  those  given 
by  equations  (1)  and  (2). 
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The  various  lift  coefficients  are  related  by  equations  (1)  and  (2),  and 
typical  operation  of  the  program  Is  as  follows.  Setting  the  airfoil 
angle  of  attack  sets  Cl,  and  setting  AC|_  specifies  Cl-j-  from  equa¬ 
tion  (1)  or  vice  versa  -  Independent  of  the  position  of  the  vortex  so  long 
as  It  lies  outside  of  the  airfoil.  Moving  the  vortex  will  change  Cls 
and  but  they  will  always  satisfy  equation  (2).  The  actual  relation 
of  these  lift  coefficients  Is  derived  from  a  careful  application  of 
Blaslus1  theorem  to  the  two-body  problem,  and  this  Is  thoroughly  described 
In  reference  4. 

It  should  also  be  noted  that  there  Is  no  nejt  drag  force  on  the  airfoil 
plus  vortex  system  as  far  as  potential  flow  Is  concerned,  and  Cqa  repre¬ 
sents  a  mutual  repelling  force  between  the  vortex  and  the  airfoil.  It 
Is  this  force  which  causes  leading  edge  slats  to  extend  automatically 
when  released  from  their  retracted  position. 

A  sample  result  Is  shown  In  figures  7a  and  7b.  Figure  7a  shows  the  result  of 
an  arbitrary  specification  of  the  vortex  location  and  strength.  It  can  be  seen 
that  this  results  In  a  severe  and  unwanted  velocity  spike  near  the  leading 
edge.  Leaving  ACL  fixed  and  adjusting  the  vortex  position  to  obtain  a 
smooth  velocity  modulation  yields  the  result  shown  In  figure  7b.  It  was  found 
that  It  was  quite  easy  to  select  the  proper  vortex  location  and  strength  In 
order  to  obtain  a  smooth  velocity  modulation.  This  represents  one  of  the 
major  features  and  justifications  for  employing  the  computer  graphics  capa¬ 
bility. 

The  following  section  describes  the  application  of  the  point  vortex  model 
to  some  representative  problems  of  multi-element  airfoil  analysis  and  design. 

6,2  Point  Vortex  Model;  Results  and  Conclusions 

The  pole  airfoil  plus  point  vortex  program  has  been  used  to  examine  a  limited 
set  of  leading  edge  design  and  analysis  problems.  In  addition  to  the  techni¬ 
cal  Information  obtained,  this  study  has  demonstrated  the  value  of  the  on¬ 
line  computer  graphics  capability  of  the  program  for  providing  the  engineer 
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with  a  basic  Intuition  regarding  leading  edge  slat  performance.  This  same 
Intuition  could  be  obtained  without  using  computer  graphics,  and  Instead 
running  a  variety  of  cases  and  then  reducing  the  data  to  provide  simi  lar  re¬ 
sults.  However,  the  efficiency  and  vividness  provided  by  the  graphics  capa¬ 
bility  appears  to  justify  Its  use. 

The  example  problems  which  follow  are  Intended  to  demonstrate  the  capability 
of  the  pole  airfoil  plus  point  vortex  program  and  evaluate  the  point  vortex 
as  a  theoretical  model  for  a  leading  edge  slat.  They  should  not  be  regarded 
as  a  comprehensive  study  of  all  of  the  possible  applications  of  the  program. 
Many  potentially  Interesting  problems  have  been  left  unexplored  -  due  pri¬ 
marily  to  the  limited  time  available  for  this  study. 

A  nine  percent  thick  cambered  Joukowskl  airfoil  (see  figure  6)  was  used  as 
the  main  airfoil  In  most  of  the  example  problems  which  follow.  This  profile 
was  felt  to  be  a  representative  geometry  for  leading  edge  slat  applications. 
Moderate  variations  In  airfoil  camber  and  thickness  failed  to  produce  any 
significant  qualitative  changes  In  the  results.  One  possible  exception  was 
that  If  the  airfoil  was  very  thick  (~20X),  the  effectiveness  of  the  point 
vortex  as  a  slat  was  not  very  good  In  that  It  was  necessary  to  locate  a 
strong  vortex  at  a  large  distance  from  the  airfoil  In  order  to  obtain  a  satis 
factory  modulation.  On  the  other  hand,  leading  edge  slats  are  not  typically 
employed  on  very  thick  airfoils  because  they  are  not  needed.  It  Is  likely, 
however,  that  much  useful  Information  could  be  gained  from  a  closer  examina¬ 
tion  of  the  effects  of  main  airfoil  geometry,  but  unfortunately  time  did  not 
permit  this  In  the  present  study. 

Example  1. 

It  Is  desired  to  obtain  a  lift  coefficient  of  3,0  from  a  given  airfoil 
(in  this  case  the  9%  thick  cambered  Joukowskl  airfoil)  with  and  without 
a  leading  edge  slat,  and  study  the  possible  justification  for  using  the 
slat.  This  example  will  also  show  how  two  airfoils  may  be  combined  to 
develop  a  total  of  3.0. 
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First, the  Joukowski  airfoil  was  generated  and  its  angle  of  attack  adjusted 
to  provide  =  3.0  (  a  =  20.8°)  with  no  vortex  present  (see  figure  8a). 
Next,  a  vortex  of  very  low  strength  (  AC^  =  0.03)  was  located  near  the 
leading  edge  to  provide  the  "best"  possible  suppression  of  the  pressure 
peak.  Due  to  the  presence  of  the  vortex  the  total  lift  coefficient  becomes 

C.  =  C,  +  AC.  =  3.0  +  0.03  =  3.03 
lT  l  L 

and  therefore  the  angle  of  attack  was  reduced  to  obtain  Cl-j-  =  3.0,  and 
the  result  is  shown  in  figure  8b.  The  vortex  strength  was  then  increased 
to  ACl  =  0.08  and  the  above  process  repeated,  each  time  adjusting  the 
vortex  position  to  obtain  the  best  possible  modulation  of  the  airfoil's 
pressure  peak  -  a  qualitative  judgement.  Sample  results  are  shown  in 
figures  8b,  8c,  and  8d.  In  each  case  the  total  lift  coefficient  is  the 
same,  =  3.0.  The  variations  of  the  force  coefficients  0^.,  C|_^, 

Cj_s,  and  Cj_  together  with  the  pressure  peak  Cpm1.n  of  the  modulated 
distribution  are  plotted  versus  a  in  figure  9. 

The  range  of  applicability  of  these  results  to  conventional  slat-airfoil 
performance  is  limited  by  the  following  consideration.  To  a  rather  crude 
approximation,  the  total  lift  of  a  typical  slat  plus  airfoil  combination 
is  usually  divided  between  the  two  elements  ’‘n  proportion  to  their  respec¬ 
tive  chords.  Slat  chords  are  rarely  as  much  as  20%  of  the  main  airfoil 
chord,  and  a  more  common  value  is  on  the  order  of  10%.  Referring  to 
figure  9,  20%  of  Ci  =3.0  gives  Ci  =0.6  and  Cl  =2.4  which 
suggests  the  limit  of  applicability  of  these  results  to  conventional  slat 
designs. 

An  interesting  cross-plot  can  be  obtained  form  the  results  of  figure  9. 
Assuming  the  the  slat  chord  is  equivalent  to  the  percentage  Cl^/Cl^., 
the  data  in  figure  9  may  be  used  to  generate  a  plot  of  Cp  .  versus 
slat  chord  required  to  provide  it,  i.e.,  from  figure  9  a  slat  chord  of 
10%  implies  Cl^  =  0.10 xCi.j  =  0.3  which  corresponds  to  Cpmin  =  -9  at 
a=  20.5°.  The  result  is  shown  in  figure  10.  This  should  not  be 
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regarded  as  a  precise  analytic  fact,  but  rather  as  a  qualitative  Indica¬ 
tion  of  the  relation  between  Cpm^n  and  the  slat  chord  required  to  obtain 
It.  It  can  be  seen  that  a  slat  chord  on  the  order  of  10%  of  the  main  air¬ 
foil  chord  provides  a  very  substantial  reduction  In  Cpm^n.  Increasing 
the  slat  chord  beyond  about  15%  produces  a  relatively  small  reduction  In 

Cpm1n* 

Example  2. 

Here  It  Is  desired  to  Increase  while  holding  Cpm_jn  on  the  airfoil 
to  a  specified  value  which  might  be  set  by  separation  criteria.  Again,  a 
9%  thick  cambered  Joukowskl  airfoil  was  used.  First,  the  airfoil  angle 
of  attack  was  set  so  that  Cpm1n  reached  the  specified  value  with  no 
vortex  present.  The  angle  of  attack  was  then  Increased  by  a  few  degrees 
and  the  vortex  strength  and  position  adjusted  so  that  Cpm^n  was  reduced 
to  the  specified  value  using  the  minimum  possible  vortex  strength.  This 
process  was  then  repeated  until  the  airfoil  angle  of  attack  reached  45° 
which  Is  somewhat  beyond  the  range  of  practical  operation.  Sample  results 
are  shown  In  figures  11a,  11b,  and  lie.  The  resulting  variations  of  the 
various  lift  coefficients  are  shown  In  figures  12  and  13  for  specified 
values  of  Cpm1n  of  -8  and  -12  respectively.  Also  shown  Is  the  Implied 
slat  chord  based  on  the  ratio  Cls/CLt  as  In  Example  1.  Using  this  data, 
a  cross  plot  of  Clj  versus  the  Implied  slat  chord  required  to  obtain  It 
was  made  and  Is  shown  In  figure  14  for  both  Cpm^n  *  -8  and  Cpm^n  »  -12. 

One  of  the  cruder  forms  of  boundary  layer  separation  criteria  says  that  the 
existence  of  attached  flow  Is  dependent  on  some  maximum  (numerical)  value 
of  Cpm1n  wh1ch  11es  somewhere  between  -8  and  -12.  A  pressure  peak  which  Is 
higher  Is  likely  to  cause  separation  -  possibly  of  the  entire  upper  surface 
of  the  airfoil  since  the  peak  Is  located  near  the  nose.  Therefore,  figures 
12  and  13  may  be  regarded  as  demonstrating  the  ability  of  a  slat  to  Increase 
the  CLmax  of  an  airfoil.  Alternatively,  these  results  demonstrate  the  abil¬ 
ity  of  a  slat  to  Increase  the  angle  of  attack  range  of  an  airfoil.  Figure 
14  shows  that  the  slat  chord  required  varies  effectively  linearly  with  the 
Increase  In  Clt,  and  figures  12  and  13  show  that  the  slat  chord  varies  al¬ 
most  linearly  with  a. 
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Example  3. 

The  prescription  of  the  strength  and  position  of  the  vortex  to  obtain  a 
"proper"  or  "best  possible"  modulation  on  the  airfoil  Is,  of  course,  to 
some  extent  subjective.  However,  It  was  found  that  If  these  parameters 
were  set  at  the  highest  angle  of  attack  at  which  the  system  was  to  be 
operated,  when  the  angle  of  attack  was  reduced  It  was  only  necessary  to 
reduce  ACL  and  not  change  the  location  of  the  vor  i  In  order*  to  obtain 
proper  modulation  at  the  lower  angles.  Conversely,  1*  the  vortex  location 
was  set  at  some  Intermediate  angle  of  attack  and  proper  modulation  obtained 
there.  It  was  Invariably  necessary  to  adjust  both  aC|_  and  the  vortex  posl 
tlon  as  the  angle  of  attack  was  Increased  in  order  to  provide  a  smooth  modu 
latlon. 

This  phenomenon  Is  Illustrated  by  using  the  9%  thick  Joukowskl  airfoil  and 
requiring  a  very  strong  pressure  modulation.  The  airfoil  was  set  at 
25°  and  ACl  and  the  vortex  position  (r,  0)  were  adjusted  to  ob¬ 
tain  a  smooth  modulation  as  shown  in  figure  15a. 

Reducing  a  to  20°  but  leaving  ACl*  r,  and  6  fixed  yielded  the  spiked 
pressure  distribution  shown  In  figure  15b.  However,  by  simply  reducing 
ACl  the  smooth  distribution  shown  in  figure  15c  was  obtained.  This 
process  was  continued  down  to  «  =  -6°,  where  at  each  interval  ACl  was 
reduced  just  enough  to  obtain  a  smooth  modulation  and  the  vortex  location 
remained  fixed  at  the  position  set  at  a  =  25°.  Some  sample  modulation 
distributions  are  shown  In  figures  16a,  16b,  and  16c  for  a  =  10°,  5°, 
and  0°,  respectively.  The  resulting  variations  In  the  lift  coefficients 
versus  a  are  shown  In  figure  17. 

Figure  17  shows  an  almost  linear  variation  of  Clt»  Cl  ,  and  CLs  with 
respect  to  a.  Since  ACl  was  se*  solely  on  the  basis  of  the  appearance 
of  the  resulting  modulated  pressure  distribution  this  might  not  be  an 
expected  result.  However,  If  the  free  stream  velocity  vector  is  con¬ 
sidered  as  being  made  up  of  two  components;  V  parallel  to  the  chord, 
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and  normal  to  the  chord,  a  simple  explanation  can  be  offered.  For 
small  a,  V«N  Is  nearly  directly  proportional  to  a,  and  It  Is  the  so- 
called  "cross-flow  component  of  velocity"  °N  which  causes  the  leading 
edge  pressure  peak  on  the  airfoil.  Therefore  since  \kN  varies  approxi¬ 
mately  linearly  with  a, and  the  Intensity  of  the  pressure  peak  does  like¬ 
wise,  It  Is  not  surprising  that  the  strength  of  the  vortex  required  to 
suppress  the  peak  also  varies  linearly  with  a. 

The  values  of  C|_s  which  result  are  somewhat  high  for  conventional  slat 
applications.  This  Is  a  consequence  of  the  relatively  strong  modulation 
asked  for  Initially  which  was  used  to  Illustrate  the  principles  de¬ 
scribed  above  as  vividly  as  possible.  Also,  the  unmodulated  pressure 
distributions  at  the  lower  angles  of  attack  do  not  really  require  any 
modulation  to  avoid  separation.  However,  as  stated  earlier,  ACL  was 
reduced  at  each  Interval  (value  of  a)  just  enough  to  remove  any  undesire- 
able  pressure  spikes  as  shown  In  figure  15,  and  therefore  some  modulation 
remained  -  In  fact,  a  negative  AC^  was  used  *t  the  lower  angles  of 
attack. 

Example  4. 

The  purpose  of  this  example  was  to  test  the  validity  of  approximating  a 
slat  by  a  single  point  vortex.  The  geometry  used  was  a  10%  thick  syirmetrl 
cal  Joukowskl  main  airfoil  with  a  cambered  Joukowskl  airfoil  whose  chord 
was  10%  of  the  main  airfoil's  for  the  slat.  The  slat  was  located  In  a 
typical  position  with  respect  to  the  main  airfoil  as  shown  In  figure  18. 
The  combination  was  then  input  Into  the  Douglas  Neumann  potential  flow 
program  at  an  angle  of  attack  of  15°  which  provided  the  aerodynamic  coef¬ 
ficients  of  the  system  and  Its  elements  together  with  the  pressure  dis¬ 
tributions  on  the  elements. 

Next,  using  the  pole  airfoil  plus  point  vortex  program,  a  point  vortex 
was  located  at  the  slat  quarter  chord  point  and  aCl  was  set  so  that  the 
total  load  on  the  vortex  (Cls2  +  Cd^2)1/2  was  equal  to  the  total  load  on 
the  slat  as  calculated  from  the  Douglas  Neumann  program  results.  The 
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resulting  modulation  distributions  from  both  the  point  vortex  and  the 
slat  are  shown  In  figure  18,  along  with  the  values  for  Clt  and  C^. 

The  slat  provides  slightly  more  modulation  than  the  vortex  and  thus  the 
total  lift  of  the  slat  plus  airfoil  system  (Clt  =  1.79)  Is  less  than  that 
point  vortex  plus  airfoil  system  (C^  =  1.84)  for  the  same  loading  of 
the  slat  and  point  vortex. 

Since  the  slat  Is  positioned  less  than  Its  chord  length  away  from  the 
nose  of  the  airfoil,  the  above  result  Is  not  surprising.  It  Is  expected 
that  the  distributed  vortlclty  along  the  slat  chord  would  provide  more 
effective  modulation  than  the  concentrated  vortlclty  of  the  point  vortex. 
However,  these  results  do  demonstrate  that  using  a  point  vortex  to  repre¬ 
sent  a  slat  as  a  first  order  theoretical  model  Is  a  reasonable  approach. 

Example  5. 

The  point  vortex  can  also  be  used  to  simulate  a  slotted  flap  at  the  air¬ 
foil  trailing  edge.  An  example  of  this  Is  shown  In  figure  19  where  the 
vortex  Is  located  just  below  the  trailing  edge.  The  circulation  about 
the  flap  (point  vortex  here)  causes  an  acceleration  of  the  flow  near  the 
trailing  edge  on  the  main  airfoil  which  reduces  the  pressure  recovery 
requirement  at  the  trailing  edge.  However,  this  also  Increases  the  vel¬ 
ocity  peak  at  the  leading  edge  as  can  be  seen  In  figure  19  which  calls 
for  a  leading  edge  slat  (or.  In  the  present  context,  a  second  point  vor¬ 
tex).  At  this  time,  the  pole  airfoil  plus  point  vortex  program  does  not 
have  the  capability  for  generating  two  independent  point  vortices,  but 
It  appears  that  this  would  be  a  logical  next  step  In  developing  the  pro¬ 
gram. 
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7.0  DISTRIBUTED  SINGULARITY  MODEL 


7.1  General  Discussion 


As  mentioned  In  the  Introduction,  a  next  level  of  sophistication  beyond  the 
point  vortex  model  appears  to  be  the  representation  of  the  slat  by  a  contin¬ 
uous  distribution  of  singularities  along  the  slat  chord.  The  use  of  thin 
airfoil  theory  for  describing  the  slat  and  Its  effect  on  the  main  airfoil  Is 
particularly  appealing  from  the  standpoint  that  linearized  theory  allows  the 
separation  of  camber  and  thickness,  and  moreover,  the  theory  Is  quite  simple 
and  straightforward  when  compared  with  the  exact  solution  formulation.  Since, 
as  discussed  earlier,  the  exact  goals  of  slat  analysis  and  design  theory 
are  not  completely  clear.  It  Is  felt  that  a  linearized  analysis  should  be 
conducted  In  order  to  provide  proper  guidance  for  an  exact  solution  study. 

It  may  be  recalled  that  much  of  the  fundamental  Information  regarding  the 
performance  of  single  element  airfoils  comes  directly  from  thin  airfoil 
theory. 

The  basic  "new"  analytical  problem  of  using  thin  airfoil  theory  to  represent 
the  slat  Is  the  determination  of  an  appropriate  system  of  Imaging  the  slat 
singularity  distributions  so  that  the  presence  of  the  main  airfoil  Is  pro¬ 
perly  accounted  for.  At  first.  It  appeared  that  an  extension  of  the  point 
vortex  model  to  the  case  of  distributed  vortlclty  along  a  slat  chord  line 
might  be  the  best  approach  to  use.  However,  examination  of  the  method  devel¬ 
oped  by  O'Pray  In  reference  4  suggested  that  continued  development  and  eval¬ 
uation  of  his  work  would  probably  be  the  most  fruitful  course  to  follow. 

Since  the  available  documentation  of  O'Pray's  method  was  rather  limited.  It 
was  first  necessary  to  redevelop  his  entire  analysis  here  at  Douglas.  Also, 
O'Pray  obtained  his  sample  solution  using  an  on-line  computer  facility  which 
allowed  him  to  directly  Interact  with  the  calculation  procedure  and  adjust 
various  Inputs  and  parameters  until  a  satisfactory  result  was  obtained.  There 
fore  It  was  decided  to  develop  a  new  computer  program  which  would  minimize 
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if  not  eliminate  any  requirements  for  "man-in-the-loop"  interaction.  It  is 
stressed  that  the  above  comments  should  in  no  way  be  interpreted  as  a  criti¬ 
cism  of  O'Pray's  work.  Instead,  his  synthesis  and  formulation  of  the  lead¬ 
ing  edge  slat  design  problem  was  regarded  as  sufficiently  impressive  that 
further  development  was  definitely  in  order.  The  discussion  which  follows 
is  intended  as  a  detailed  description  but  not  a  rigorous  development  of  O'Pray's 
method. 

The  problem  to  be  solved  is  that  of  designing  a  leading  edge  slat  which  pro¬ 
vides  a  specified  pressure  distribution  modulation  on  the  nose  region  of  the 
main  airfoil.  As  stated,  this  problem  is  somewhat  underspecified  and  it  is 
not  likely  to  possess  a  unique  solution.  Therefore  some  additional  conditions 
may  be  imposed.  The  use  of  thin  airfoil  theory  implies  representing  the  slat 
by  an  infinite  set  of  singularity  distributions  whose  coefficients  are  to  be 
set  on  the  basis  of  the  boundary  conditions  which  are  to  be  satisfied.  The 
pressure  field  imposed  on  the  main  airfoil  by  these  singularity  distributions 
is  linear  with  respect  to  the  magnitude  of  their  coefficients  and  nonlinear 
with  respect  to  the  location  of  the  slat  chord.  Thus  it  seems  reasonable, 
at  least  to  begin  with,  to  specify  the  slat  chord  position,  and  then  attempt 
to  determine  the  proper  values  for  the  mode  coefficients.  The  choice  of  the 
coefficients  involves  two  basic  requirements:  (1)  they  must  provide  the 
desired  pressure  distribution  modulation  on  the  main  airfoil,  and  (2)  they 
must  correspond  to  a  real  and  useable  geometry  for  the  slat  itself.  (There 
is  no  a  priori  guarantee  that  an  arbitrarily  chosen  set  of  mode  coefficients 
might  not  imply  a  slat  of  negative  thickness  or  one  whose  trailing  edge  does 
not  close.) 

One  of  the  major  features  of  O'Pray's  method  lies  in  the  choice  of  the  various 
conformal  mappings  employed.  Instead  of  working  with  an  arbitrary  main  air¬ 
foil  shape,  an  ellipse  of  arbitrary  thickness  is  used.  It  is  observed  that  for 
most  conventional  airfoils  without  extreme  leading  edge  camber  or  nose  droop, 
the  nose  region  itself  is  very  nearly  elliptical,  and  therefore  the  nose  flow 
about  an  ellipse  closely  models  the  typical  flow  environment  of  a  leading  edge 
slat.  Thus  an  ellipse  with  its  circulation  set  to  maintain  the  Kutta  condition 
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at  Its  trailing  edge  becomes  a  rather  Ideal  model  for  a  main  airfoil  to  study 
slat  performance  since  Its  mapping  to  a  circle  Is  straightforward. 

7.2  Description  of  Analysis. 


A.  Mapping  functions. 

Several  consecutive  mappings  are  used  In  the  analysis  beginning  with  the 
mapping  from  the  ellipse  (  Z  (x-j,  yj)-plane,  figure  20a)  to  the  circle 
(z  (x2»  y2)-plane,  figure  20b)  where  the  mapping  function  Is  given  by 


Z  ■  z  +  C7/4Z 
2  *  V2[z  +  (Z2  -  c7)‘/2] 


(3) 


where  Cy  a  a2  -  b2  and  a  and  b  are  the  major  and  minor  axes  of  the 
ellipse.  The  radius  c3  of  the  corresponding  circle  In  the  z-plane  Is  ob¬ 
tained  from  c3  a  (a  +  b)/2.  Next,  the  circle  Is  opened  at  Its  trailing 
edge  and  Is  mapped  to  the  real  axis  of  the  half-plane  (W(x3,  y3)-plane, 
figure  20c)  where  y3  >  0  corresponds  to  that  region  outside  the  circle. 
The  point  at  Infinity  In  the  z-plane  maps  to  y3  ■  1  In  the  W-plane,  and 
the  leading  edge  maps  to  the  origin  while  the  leading  edge  stagnation  point 
maps  to  a  small  value  of  -  X3  on  the  negative  real  axis  (If  the  angle  of 
attack  of  the  ellipse  Is  zero,  the  stagnation  point  and  the  leading  edge 
become  coincident).  The  mapping  from  the  circle-plane  to  the  half-plane  Is 


W 

z 


(4) 


The  genr.  1e?  of  these  mappings  are  shown  In  figure  20,  where  the  modulation 
re"-'’ 1'  Indicated  by  the  heavy  line.  Also  shown  In  figure  20c  are  the  stream' 
Hn.s  0;  ,jie  flow  near  the  leading  edge.  It  Is  observed  that  In  the  near 
vicinity  of  the  modulation  region  the  curvature  of  the  streamlines  Is  rather 
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moderate,  and  therefore  a  straight  line  slat  chord  may  be  placed  tangent  to 
one  of  the  streamlines  as  shown  In  figure  20c.  The  selection  of  a  particular 
slat  chord  position  Is  based  on  the  location,  extent,  and  character  of  the 
desired  modulation  distribution.  This  will  be  discussed  In  greater  detail 
later. 

Once  the  slat  chord  position  Is  set,  a  simple  mapping  to  the  slat-centered 
half-plane  (s  (s,t)  -  plane)  Is  made,  where  the  origin >of  the  S-plane 
lies  at  the  slat  midchord  and  the  real  axis  parallels  the  slat  chord  as  shown 
In  Figure  20c.  This  mapping  Is  merely  a  translation  plus  rotation  plus  dila¬ 
tion  and  is  given  by 


S 


(5) 


where  f2  and  f-j  are  the  coordinates  of  the  slat  midchord,  a%  Is  the  slat 
angle  with  respect  to  the  x3-axis,  and  c$  Is  the  slat  chord  length  -  all 
measured  in  the  W  -plane.  Since  It  has  been  assumed  that  the  slat  chord  will 
be  set  tangent  to  the  local  streamline  which  determines  «s,  only  the  three 
parameters  fj,  f2,  and  c$  are  needed  to  specify  the  slat.  In  the  S-plane, 
the  slat  chord  runs  from  s  =  -2  (the  slat  leading  edge)  to  s  =  +2  (the  slat 
trailing  edge)*.  This  completes  the  mappings  used  in  the  analysis. 

The  main  appeal  of  the  above  set  of  mappings  lies  In  the  fact  that  when  a 
singularity  distribution  Is  defined  on  the  slat  chord,  an  Identical  singular¬ 
ity  distribution  may  be  defined  on  an  Image  slat  chord  located  as  shown  in 
figure  20c.  This  Is  the  same  procedure  used  for  analyzing  a  conventional  air¬ 
foil  In  ground  effect  where  the  ground  surface  becomes  the  axis  of  symmetry 
of  the  potential  flow  field.  The  fact  that  the  main  flow  field  in  the  W  - 
plane  Is  not  a  uniform  stream  does  not  affect  this  model.  In  the  analysis 
which  follows,  the  X3  -  axis  of  the  W-plane  (which  is  actually  the  surface 


*A  slat  chord  of  4  is  convenient  in  defining  the  singularity  distributions. 
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of  the  ellipse),  will  often  be  referred  to  as  the  "ground  plane"  In  analogy 
to  the  above. 

B.  Singularity  distributions. 

The  singularity  distribution  used  to  represent  the  velocity  field  Induced  by 
the  slat  comes  from  classical  thin  airfoil  theory  and  has  the  form 

q  (()  "  {+T  +  »  cn  "  an  +  1  bn  <6) 

where  q  ■  u  -  iv  Is  the  complex  velocity  and  (  Is  related  to  the  slat-plane 
variable  S  by  the  transformation 

S  -  {+  1/C,  £  -  S-+JS2  -  41  /2  (7) 


It  should  be  noted  that  this  Is  not  a  mapping  function,  but  rather  an  Inde¬ 
pendent  variable  transformation,  l.e., 

q  (S)  •  q  [{  (S)l 

and  therefore  In  terms  of  S,  equation  (6)  becomes 


q  (S) 


[Lzisi-ti 

2 


(8) 


The  use  of  (  as  an  Independent  variable  Is  merely  for  convenience.  Expanding 
the  coefficients  of  equation  (6)  yields 


q  (0 


ao 

FT 


A 

1  c" 


1  bft  jo  bn 

+  I+T+  jif 
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where  the  an's  represent  slat  thickness,  and  the  bn's  represent  slat 
camber.  In  order  for  the  corresponding  slat  to  close  at  Its  trailing  edge 
the  requirement 


-  0 


must  be  satisfied,  which  requires  aQ  +  a^  ■  0.  For  a  finite  nose  radius 
aQ  >  0,  and  therefore  a-|  Is  given  by  a-j  *  -  aQ. 


In  the  case  of  the  slat  design  problem  considered  here,  the  flow  field  Im¬ 
posed  on  the  main  airfoil  by  the  slat  will  not  be  significantly  affected  by 
the  higher  order  terms  of  the  singularity  distribution.  Moreover,  Inclusion 
of  a  large  number  of  terms  In  the  least  squares  matching  analysis  (the  discus¬ 
sion  of  which  follows  shortly)  would  tend  to  create  numerical  difficulties  since 
the  details  carried  by  the  higher  order  terms  are  not  readily  Identifiable  when 
one  moves  off  the  slat  surface  Itself.  Therefore  the  singularity  distribu¬ 
tion  Is  defined  as  consisting  of  the  first  four  thickness  modes  plus  the  first 
four  camber  modes.  Since  the  first  two  thickness  modes  are  coupled  by  the 
closure  condition  (a^  *  -aQ),  the  total  number  of  Independent  modes  Is 
actually  seven. 


Using  the  form  of  equation  (8)  and  redefining  the 
the  seven  singularity  modes  are  defined  as: 


an's  and  bn's 


camber  modes 


as 


Bj's, 


(9) 
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The  corresponding  slat  geometry  mode  shapes  are  sketched  In  figure  21,  and 
the  velocity  distribution  of  equation  (6)  becomes 


q(S) 


Bj  qj(S) 


In  terms  of  S.  The  complex  velocity  In  the  half  plane  q  (W)  Is  given  by 


q(W)  ■  q(S)  ■  q(S)£-e^*s 


(10) 


where  dS/dW  Is  obtained  from  equation  (5).  Now  the  real  and  Imaginary  parts 
of  q(W)  on  W  =  X3  (y3  =  0)  represent  the  parallel  and  normal  slat  singular¬ 
ity  Induced  velocities  on  the  ground  plane  (l.e.  the  ellipse  surface  In  the 
W-plane).  Defining 


qj(S)  *  Uj(S)  —  1  Vj(S) 

qj(w)  a  uj (w)  -  1  Vj(w) 

and  using  equation  (10)  yields 

UI (W)  =  [uj(S)  cos  «s  +  Vj(S)  sin  agJ 

VI (W)  *  j^Uj(S)  sin  a$  -  Vj(S)  cos  a 


s] 


(11) 
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as  the  components  of  the  velocity  Induced  In  the  W-plane  (half-plane)  by 
each  of  the  I  singularity  modes,  where  the  Uj(S)  and  Vj (S )  are  known 
from  the  expressions  given  by  equations  (9). 


Since  an  Image  slat  is  located  In  the  lower  half  of  the  W-plane,  the  x3 
axis  becomes  an  axis  of  symmetry,  and  on  the  X3  -  axis  the  total  velocity 
Induced  by  the  slat  plus  Image  slat  singularity  distributions  Is  simply 


u<*3> 


B,  u, 


(12) 


where  the  uj(W)  are  defined  by  equations  (11),  and  are  functions  of  the 
slat  chord  and  Its  position  In  the  W-plane  only.  The  factor  2  Is  due  to 
the  presence  of  the  Image  slat. 


C.  Determination  of  the  mode  coefficients. 


The  problem  remains  to  determine  the  coefficients  Bj  of  the  singularity 
distribution  modes  so  that  the  Induced  velocity  u(x3)  provides  the  desired 
modulation  on  the  ground  plane  (l.e.  on  the  nose  region  of  the  ellipse).  Typi¬ 
cally  the  desired  velocity  modulation  will  be  specified  by  a  set  of  values  at 
a  finite  number  of  discrete  points  along  the  ellipse  nose  surface.  This 
specified  velocity  Is  readily  mapped  to  the  half-plane  where  It  Is  defined  as 


w (x3i )  =  wi*  1=1,2,  — ,  L  (number  of  matching 

stations) 


where  the  x^  are  the  matching  stations.  Next,  comparing  the  desired 
velocity  with  that  Induced  by  the  singularity  distributions  yields  the 
L-equatlons 


a. 


w(x3i)  “  2  S  BJ  UJ  (X31  ^  *  *  =  L 

J=  1 


or 


"1-2 


,  1  -  1,  L 
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The  terms  u^  are  called  "Influence  coefficients"  and  are  defined  as  the 
velocity  on  the  ground  plane  at  the  matching  station  x3i  Induced  by  the  Jth 
singularity  mode  with  a  unit  mode  coefficient. 

The  values  of  the  mode  coefficients  Bj  are  determined  using  a  standard  least 
squares  matching  procedure.  Letting 

£  3  is,2 

1=1  1 

and  differentiating  <f>  with  respect  to  the  Bj  yields 


as  seven  simultaneous  equations  to  be  solved  for  the  seven  mode  coefficients 
Bj.  In  matrix  notation,  equations  (13)  become 


u  (w  -  uT  B) 


where  the  matrices  of  equation  (14)  are  defined  by 


u 

=  1 

1 

uji  1 ' 

J  x  1  rectangular  matrix 

w 

,  1 

I 

*1  !  * 

1  column  matrix 

B 

■  ! 

bj  I  ■ 

J  column  matrix 
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Equation  (14)  can  be  rewritten  in  the  form 
uuTB  =  u  w 

or  A  B  =  C  (15) 

where  A  =  u  u\  and  C  =  u  w.  Equation  (15)  is  a  standard  form  and  may 
be  solved  for  B  since  A  and  C  are  known. 


The  resulting  mode  coefficients  Bj  provide  the  best  possible  fit  by  the 
singularity  induced  velocity  distribution  given  by  equation  (12)  to  the  spec¬ 
ified  modulation  distribution  w(xj).  In  principle,  an  arbitrarily  large 
number  of  singularity  modes  could  be  used  for  this  calculation,  however  it 
is  likely  that  numerical  instability  would  appear  in  the  matrix  inversion. 

D.  Corrections  and  adjustments  to  the  mode  coefficients. 

Before  proceeding  to  the  determination  of  the  slat  shape  corresponding  to  the 
mode  coefficients  BJt  it  is  necessary  to  consider  some  adjustments  of  the 
Bj  which  will  improve  the  accuracy  of  the  solution.  These  are  the  adjust¬ 
ment  of  the  slat  chord  angle  and  a  correction  to  the  circulation  about  the 
ellipse. 


It  is  recalled  that  the  slat  chord  is  initially  positioned  tangent  to  the 
local  streamlines  at  midchord  which  in  effect  corresponds  to  a  zero  angle 
of  attack  in  a  uniform  free  stream.  Referring  to  figure  21,  the  B-j  and  B^ 
camber  modes  serve  to  adjust  the  slat  chord  angle  by  increasing  it  when  they 
are  positive.  It  can  be  shown  that  the  change  in  angle  is  given  by 


where  umf  is 
chord  (V^  in  a 
essarily  small 


the  local  main  flow  velocity  parallel  to  the  slat  at  its  mid¬ 
uniform  free  stream).  It  turns  out  that  A a$  is  not  nec- 
and  therefore  the  initial  slat  chord  angle  should  be  adjusted 
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accordingly.  This  changes  the  values  of  the  Influence  coefficients  and  there¬ 
fore  a  repeat  of  the  least  squares  solution  Is  required  to  obtain  a  refined 
set  of  mode  coefficients  Bj.  Repeating  this  procedure  again  Is  not  justified 
on  the  basis  of  the  Inherent  accuracy  of  thin  airfoil  theory. 

Another  correction  arises  due  to  the  Image  slat  vortlclty.  When  mapped  back 
to  the  ellipse  plane*  the  Image  slat  lies  Inside  the  ellipse  and  Its  vortlclty 
due  to  the  camber  modes  Bj,  Bp,  B^,  and  B^  moves  the  rear  stagnation  point 
and  the  Kutta  condition  Is  no  longer  satisfied.  The  Image  slat  circulation 
Is  calculated  by  considering  the  problem  In  the  circle  plane.  Most  of  the 
slat  vortljclty  Is  contained  In  the  Bj  and  Bp  terms,  and  the  B^  and  B^ 
terms  can  be  neglected  within  the  accuracy  of  thin  airfoil  theory.  The  vor¬ 
tlclty  of  the  and  Bp  terms  Is  approximated  by  a  point  vortex  of  strength 
2 ttBj  at  the  slat  quarter  chord  and  a  point  vortex  of  strength  27rBp  at  the 
slat  midchord.  The  corresponding  Image  locations  and  circulation  correction 
required  to  maintain  the  Kutta  condition  Is  then  calculated  In  the  circle  plane 
using  the  same  method  as  was  employed  In  the  point  vortex  model. 

Even  though  the  ellipse  circulation  correction  amounts  to  less  than  one  percent 
of  the  total  ellipse  circulation,  It  causes  a  significant  change  In  the  ellipse 
pressure  distribution,  particularly  In  the  nose  region.  As  a  result  the  spec¬ 
ified  modulation  velocity  distribution  w(x^)  changes  and  therefore  a  new 
least  squares  solution  Is  necessary.  Due  to  the  sensitivity  of  the  nose  flow 
to  the  ellipse  circulation  this  procedure  Is  repeated.  Convergence  to  a  negli¬ 
gible  correction  occurs  after  three  or  four  cycles. 

E.  Slat  thickness  constraints. 

Thus  far  In  the  analysis  the  calculation  of  the  mode  coefficients  has  been 
based  almost  solely  on  the  requirement  that  they  provide  the  desired  velocity 
modulation  distribution  on  the  nose  of  the  ellipse.  It  now  becomes  necessary 
to  Identify  some  geometrical  considerations  -  In  particular  with  respect  to  the 
slat  thickness  distribution.  Earlier  In  the  analysis  closure  of  the  slat  was 
guaranteed  by  setting  a-j  *  -  aQ.  This  does  not,  however.  Insure  that  the 
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thickness  mode  coefficients  Bg,  Bg,  and  By  as  chosen  by  the  least  squares 
solution  won't  provide  a  slat  thickness  distribution  which  is  re-entrant.  In 
fact,  it  turns  out  that  if  the  solution  is  left  unconstrained  this  is  quite 
often  what  happens.  Such  a  result  is  readily  understood  when  it  is  recalled 
that  the  slat  is  called  on  to  provide  a  velocity  modulation  which  opposes  the 
main  flow  on  the  ellipse.  This  is  accomplished  by  the  circulation  generated  by 
the  slat,  however,  in  many  cases  positive  thickness  will  tend  to  reduce  this 
effect  by  causing  an  acceleration  of  the  flow. 

This  difficulty  is  overcome  by  simply  removing  the  thickness  mode  coefficients 
from  the  least  squares  matching  solution,  and  prescribing  them  to  provide  any 
thickness  distribution  desired.  The  only  modification  required  for  the  least 
squares  solution  is  that  the  differences  become 


8i  =  w^x3i)  “wt(x3i}  “2  S  BJ  UJ  (x3i}*  1  =  L 

<J=  I 


where  wt(Xg..)  represents  the  velocity  induced  on  the  ground  plane  at 
by  the  prescribed  slat  thickness  distribution,  and  of  course  there  are  only 
B-j ,  B2,  B^,  and  B^  (the  camber  mode  coefficients)  to  be  determined.  It  is 
convenient,  but  not  necessary,  to  specify  the  thickness  distribution  by  select¬ 
ing  values  of  Bg,  Bg,  and  By,  however  any  desired  thickness  distribution  is 
in  principle  admissible. 

For  the  examples  worked  in  this  study  a  basic  thickness  form  made  up  of  the 
Bg  and  By  thickness  modes  was  used.  It  can  be  shown  that  the  Bg  mode  coef¬ 
ficient  is  related  to  the  maximum  thickness  (of  that  mode)  by 

-  Umf  ts 

B5  -  "0755“ 

where  umf  is  the  local  main  flow  velocity  at  the  slat  midchord  and  t  is 
the  maximum  thickness  to  chord  ratio.  In  order  to  thicken  the  back  half  of 
the  slat,  a  negative  By  mode  coefficient  was  chosen  whose  magnitude  was  20% 
of  Be,  i.e., 

«  n  _  n  orv  n 


-  0.20  Be 


Bg  was  set  equal  to  zero. 
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The  apparent  necessity  for  prescribing  the  thickness  distribution  Is  not  as 
bad  as  It  may  seem  at  first.  In  design  problems,  structural  considerations 
often  dictate  thickness  distribution  requirements.  The  ability  to  specify 
the  thickness  distribution  then  allows  an  additional  degree  of  freedom  In 
the  slat  design  problem  since  the  camber  distribution  will  then  be  calculated 
to  provide  the  desired  velocity  modulation  using  the  given  thickness  distribu¬ 
tion. 

F.  Determination  of  slat  geometry. 

Once  the  camber  mode  coefficients  B-j,  B2,  B3,  and  B^  have  been  calculated. 

It  remains  to  calculate  the  corresponding  slat  shape.  Working  In  the  S-plane, 
the  standard  flow  tangency  relation  Is 

dt  v 

cTs  u 

which  upon  Integration  gives 

t(s)  ■  /  &  *  ■  /  l  ds 

+2  +2 

where  t  Is  the  slat  surface  ordinate.  The  Integration  Is  begun  at  the  slat 
trailing  edge  (s  *  +2),  and  procedes  forward  (on  the  slat)  to  the  slat  lead¬ 
ing  edge  (s  =  -2). 

This  Integral  Is  more  easily  evaluated  If  the  transformation  given  by  equation 
(7)  Is  used  which  gives  on  the  unit  circle 

s  ■  2  cos  0,  ds  »  -2  s1n0d0,  ({*  e^) 

and  hence  the  Integral  takes  the  form 

e 

t(s)  =  -  2  J  ~  (0)  sin  0d0 
o 
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This  Is  simplified  to 

0 

t(s)  ■  y*v  (0)  sin  0d0  (17) 

o 

using  the  standard  thin  airfoil  theory  approximation  that  u  Is  constant 
over  the  slat  chord. 


The  normal  velocity  distribution  along  the  slat  chord  v(0)  Is  made  up  of 
three  basic  contributions:  v$(0)  due  to  the  singularity  distributions  on 
the  slat  chord,  v^m(0)  due  to  the  Image  slat,  and  vm(0)  due  to  the  main 
flow.  More  specifically,  v$(0)  Is  the  Imaginary  part  of  q  ({)  or  q(S), 
and  using  equations  (9)  It  has  the  form 

v$(0)  *  -  B1  -  B2  cos  0  +  B3  cos  20  -  cos  30 

s1n0^-  Bg  sin  20  +  B 

where  for  the  bracketed  terms,  the  (+)  applies  to  the  upper  surface  of  the 
slat  and  the  (-)  applies  to  the  lower  surface.  When  substituted  Into  equation 
(17),  vs(0)  can  be  Integrated  In  closed  form. 

The  distribution  v^m(0)  represents  the  normal  velocity  Induced  at  the  slat 
by  the  singularity  distributions  on  the  Image  slat.  vm(0)  Is  the  normal 
velocity  distribution  along  the  slat  chord  due  to  the  fact  that  the  main  flow 
streamlines  In  the  slat-plane  are  actually  curved.  If  this  were  not  Included, 
the  calculated  slat  shape  would  be  erroneously  cambered  by  an  amount  equivalent 
to  the  streamline  curvature.  Both  v.|m(0)  and  vm(0)  must  be  Integrated 
numerically  In  equation  (17). 

Once  the  geometry  solution  t(s)  Is  obtained.  It  Is  mapped  back  to  the  el¬ 
lipse  plane  using  the  mapping  functions  of  equations  (5),  (4),  and  (3). 
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G.  Method  of  calculation. 


In  order  to  evaluate  the  Inverse  slat  design  method  just  described,  a  computer 
program  was  written  In  Fortran  IV  language  for  the  IBM  System  360  computer.  The 
required  Input  data  Includes:  the  ellipse  thickness  ratio,  the  ellipse  angle 
of  attack,  the  desired  pressure  modulation  distribution,  the  thickness  dis¬ 
tribution  of  the  slat,  the  slat  chord  length  c$,  and  the  slat  midchord  loca¬ 
tion  (f2,  f i ) .  The  program  selects  the  proper  slat  chord  angle  as  and  deter¬ 
mines  the  camber  mode  coefficients  B-j,  B  ,,  I  ( ,  and  B^  which  provide  the 
closest  possible  match  to  the  desired  pressure  modulation  distribution.  These 
calculations  Include  the  corrections  for  Image  Induced  circulation  on  the 
ellipse  and  slat  chord  angle.  The  resulting  singularity  distributions  are 
then  Integrated  to  obtain  the  slat  geometry  Including  the  necessary  adjustments 
for  the  presence  of  the  Image  slat  and  streamline  curvature  In  the  W-plane. 
Finally,  the  resulting  geometry  Is  mapped  back  to  the  ellipse  plane  to  yield 
the  actual  slat  shape. 

In  Its  present  form,  the  program  requires  that  the  desired  modulation  distri¬ 
bution  be  Input  as  a  function  of  the  Z-plane  variable  x^,  and  of  course 
the  Input  variables  c$,  fj,  and  f2  are  defined  In  terms  of  W-plane  coor¬ 
dinates.  This  was  done  for  the  sake  of  expediency  In  order  to  be  able  to 
evaluate  the  distributed  singularity  distribution  model  In  the  limited  time 
available.  As  a  consequence,  prescribing  this  Input  data  is  somewhat  more 
awkward  than  would  be  desirable  -  particularly  for  an  Individual  with  limited 
experience.  However,  once  this  data  has  been  specified,  the  calculation  of 
the  slat  geometry  Is  completely  automatic  and  requires  no  Interaction  with  a 
"man-ln-the-loop".  The  computer  time  required  for  a  complete  run  on  a  360/65 
computer  Is  less  than  one  tenth  of  a  minute,  and  therefore  Iterating  to  obtain 
a  desired  slat  geometry  Is  quite  economical. 

It  should  be  noted  that  the  program  always  yields  a  slat  shape.  However,  If 
Inappropriate  values  of  c$,  f-j ,  or  f2  are  specified,  or  an  "Impossible" 
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modulation  is  sought,  the  resulting  slat  geometry  will  likely  be  rather  pecu¬ 
liar,  and  its  modulation  distribution  may  not  completely  agree  with  the  one 
requested.  Some  sample  results  are  presented  In  the  following  section. 

7.3  Distributed  Singularity  Model:  Results  and  Conclusions 

The  O'Pray  Inverse  slat  design  method  has  been  applied  to  a  limited  set  of 
examples  using  the  computer  program  developed  In  this  study.  These  examples 
are  Intended  to  illustrate  the  basic  viability  of  the  method,  and  should  not 
be  regarded  as  a  comprehensive  display  of  all  possible  applications. 

The  main  airfoil  ellipse  used  for  all  of  the  examples  Is  that  chosen  by  O'Pray 
for  his  single  test  case:  a  14%  thick  ellipse  which  In  turn  has  a  nose  radius 
of  about  one  percent  of  Its  chord.  The  ellipse  Is  set  at  an  angle  of  attack 
of  about  17.2°  (actually  0.3  radians)  which  gives  a  lift  coefficient  of  2.12, 
and  a  leading  edge  pressure  peak  of  Cp^  =  -22.7.  A  pressure  distribution 
modulation  Is  specified  which  reduces  this  pressure  peak  to  Cpm^n  =  -9.  Fig¬ 
ure  22  shows  the  pressure  distribution  on  the  nose  region  of  the  ellipse  with¬ 
out  modulation  and  that  which  results  from  the  specified  modulation  distribu¬ 
tion.  (The  other  data  shown  In  this  figure  will  be  explained  shortly.)  This 
pressure  peak  and  required  modulation  of  it  are  quite  severe,  and  thus  should 
represent  a  good  test  for  the  method. 

The  first  problem  considered  was  similar  to  O'Pray's  test  case:  design  a  10% 
thick  slat  to  provide  the  specified  modulation.  Using  the  basic  thickness 
distribution  described  In  section  7.2E  with  ts  =  0.10,  the  slat  shown  In 
figure  23a  was  obtained.  This  slat  geometry  plus  the  ellipse  was  Input  into 
the  Douglas  Neumann  potential  flow  program  and  the  modulated  pressure  distri¬ 
bution  provided  by  the  slat  Is  plotted  In  figure  22.  The  resulting  slat  geo¬ 
metry  does  not  quite  agree  with  that  obtained  by  O'Pray,  however,  this  result 
Is  not  unexpected  since  his  calculation  routine  was  not  the  same  as  that  of 
the  Douglas  program. 
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Next,  the  slat  thickness  was  Increased  to  152  by  setting  tc  ■  0.15*  and  to 
20%  by  setting  t$  «  0.20.  The  resulting  slat  geometries  are  shown  In  figures 
23b  and  23c,  respectively,  and  the  resulting  modulated  pressure  distributions 
as  obtained  from  the  Douglas  Neumann  program  are  shown  In  figure  22. 

Finally,  the  ellipse  angle  of  attack  was  Increased  to  a  «  18.4°  (0.32  radians) 
which  gives  a  lift  coefficient  of  2.3  (Cp^  ■  -25.7),  and  the  same  modulation 
(Cpmin  ■  -9)  requested.  A  15%  thick  slat  was  called  for  and  Its  resulting 
geometry  Is  shown  In  figure  23d  with  the  Douglas  Neumann  calculation  of  Its 
corresponding  modulated  pressure  distribution  again  plotted  In  figure  22. 

The  results  shown  In  figures  22  and  23  Indicate  that  there  Is  no  unique  sol¬ 
ution  to  problem  (b)  as  described  In  the  Introduction,  and  thus  there  exists 
a  wide  variety  of  slat  geometries  capable  of  accurately  matching  a  specified 
pressure  modulation  distribution  on  the  nose  region  of  a  given  airfoil.  This 
means  that  some  of  the  slat  geometry  parameters  (e.g.  the  thickness  distri¬ 
bution)  may  be  used  to  satisfy  other  constraints.  These  Include  structural, 
mechanical,  and  aerodynamic  considerations.  For  example.  It  appears  that  It 
should  be  possible  to  constrain  the  upper  surface  shape  of  the  slat,  and  with 
a  proper  reformulation  of  the  method  obtain  solutions  to  problem  (a)  described 
in  the  Introduction.  Similarly,  It  Is  possible  that  the  method  could  be  ex¬ 
tended  to  handle  the  design  problem  where  the  upper  surface  velocity  distri¬ 
bution  would  be  constrained  to  avoid  boundary  layer  separation  on  the  slat 
Itself.  Unfortunately,  the  limited  time  available  for  this  study  did  not 
permit  such  further  Investigations. 

On  the  basis  of  the  limited  application  and  results  obtained,  the  distributed 
singularity  model  as  developed  by  O' Pray  In  reference  3  appears  to  be  a  prac¬ 
tical  slat  design  tool.  At  this  time  the  method  Is  still  somewhat  embryonic 
and  requires  some  further  refinement  before  It  can  be  considered  operational. 
This  would  Include  the  ability  to  Input  arbitrary  main  airfoil  geometries 
and  their  pressure  modulations  In  the  physical  plane. 
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Operation  of  the  present  program  is  straightforward  once  one  becomes  used  to 
working  in  the  W-plane.  The  selection  of  the  proper  values  of  cs,  f  ^ ,  and 
f2»  is  based  on  the  character  and  location  of  the  modulation  distribution 
in  the  W-plane.  A  poorly  chosen  value  for  one  of  these  parameters  will 
usually  provide  a  reasonable  match  of  the  modulation  distribution,  however, 
the  resulting  slat  geometry  may  be  unacceptable.  For  example,  a  value  of  cg 
which  is  too  small  will  result  in  a  slat  which  is  highly  cambered  at  both 
ends  (i.e.  a  large  value  of  B^)  since  this  will  tend  to  broaden  the  induced 
field.  The  value  of  f2  will  tend  to  shift  the  camber  from  one  end  to  the 
other  of  the  slat.  In  general,  it  has  been  found  that  by  examining  the  agree 
ment  of  the  modulation  distributions  together  with  the  values  of  the  mode 
coefficients,  an  acceptable  slat  design  can  be  obtained  in  less  than  five 
iterations.  Since  the  total  run  time  for  the  program  on  an  IBM  System  360/65 
computer  is  less  than  0.10  min.,  this  procedure  is  quite  economical. 
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AIRFOIL  GEOMETRY  -  SLAT  RETRACTED 


AIRFOIL  GEOMETRY  -  SLAT  EXTENDED 


Figure  2. -Slat  designed  to  retract  to  form  main  airfoil  leading  edge. 
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Figure  3. -Flow  about  a  circular  cylinder  with  circulation  rr. 
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Figure  5-Circ^piane  used  for  defining  pole  airfoil  and  initial  point  vorte* 
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GENERAL  POLE  DEFINED  AIRFOIL 
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Figure  6.-Pole a  rfoil  geometry  and  velocity  distributions  on  airfoil  alone  and 
surface  beoinni^a/thi'?*'  S  ls,the  arc  ,en9th  around  the  airfoil 
at rfoi  1  perimeter  equals  iST  SUrfaCe  trai,1"s  ed9e-  The  total 
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SLAT  CHORD  IN 
%  MAIN  AIRFOIL  CHORD 

for  Example  1.  Cpm^n  versus  implied  slat  chord  required  to 
t.  CLt*  3.0. 
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Results  for  Example  2.  Vortex  strength  and  position  adjusted 
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Figure  15. -Sample  results  for  Example  3.  (a)  vortex  strength  and  position  set  to 
obtain  proper  modulation  at  a  -  25°,  (b)  alpha  reduced  to  20°  vortex 
same  as  (a),  and  (c)  vortex  strength  adjusted  to  obtain  proper  modulation 
with  vortex  position  the  same  as  (a). 
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Figure  16. -Sample  results  for  Example  3.  Vortex  position  fixed  at  «  =25°. 
Vortex  strength  adjusted  to  obtain  proper  modulation  as  a  is 
reduced:  (a)  a  =  10°,  (b)  a  =  5°,  and  (c)  a  =  0°. 
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GENERAL  POLE  DEFINED  AIRFOIL  -  VELOCITY  PLANE 

CASB  JSN  ALPHA  a  IS. 00  CL  *  I.TSSI 


CEL  CL  >  0.0000  CO  AIAF  s  0.0460 


Figure  18. -Results  for  Example  4.  Comparison  of  point  vortex  with  a  real  slat 
using  the  Douglas  Neumann  potential  flow  program. 
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Figure  19. -Results  for  Example  5.  Point  vortex  used  to  simulate  a  slotted  flap 
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Figure  20. -Domains  used  in  the  distributed  singularity  analysis;  (a)  ellipse- 
plane,  (b)  circle-plane,  and  (c)  half-plane.  (Reference  3). 


e, 


b2 


B 


3 


b4 


Figure  21  .-Geometric 
( Reference 


MDC  J5195 


CPmi„s-25  7 


-20  h 


,F>£  C  ISTRIBUTION 


CALCULATED  BY 

_4  |_  DOUGLAS  NEUMANN  PROGRAM 
O  10%  THICK  SLAT  \ 
a  15%  THICK  SLAT  i  <*=17.2° 
A  20%  THICK  SLATj 
O  15%  THICK  SLAT,  a  =  18. 4° 


_L 


-2.0 

NOSE  OF 
ELLIPSE 


-1.96 


-1.92 


_l _ 

-1.88 


I 

-1.84 

ELLIPSE 

CH0RD=4.0 


Figure  22-->;]]JPse  "ose  pressure  modulation:  desired  modulated  pressure 
distribution  and  that  provided  by  the  slats  of  figure  23. 
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designed  to  provide  the  same  modulation 


